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Coulomb matrix elements, and by inclusion of higher subbands along with several Landau levels in 
the Hilbert space. Detailed analysis of the accuracy and convergence of the exact diagonalization 
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I. INTRODUCTION 

Photoluminescence (PL) is a powerful method in the 
experimental studies of two-dimensional (2D) electron 
systems in high magnetic fieldsi 2 * 3 - In a typical exper- 
iment, an additional electron-hole (e— h) pair is intro- 
duced into the electron system through absorption of 
a photon. This causes energy relaxation via Coulomb 
scattering, followed by optical recombination of the well- 
defined, quantized eigenstates involving an e-h pair cou- 
pled in different ways to the surrounding electrons. 

This coupling is a response of the 2D electron system to 
the perturbation associated with the photo-injection of a 
hole (particle of the opposite charge). Depending on the 
nature of this response and on the number of electrons in- 
volved, the PL spectrum either reveals information about 
dynamics of an illuminated system^ or "merely" probes 
the few-body excitonic complexes which form under given 
conditions (quantum well geometry, magnetic field, etc.) 
but are weakly affected by the surrounding electrons. 

The latter scenario occurs at sufficiently high magnetic 
fields B and low electron concentrations g, correspond- 
ing to small Landau level (LL) filling factors v — 2irg\ 2 
(with A = y/hc/eB being the magnetic length). The PL 
spectra at v <§C 1 are usually dominated by recombination 
of three-body states called trions (X~ = 2e + h), consist- 
ing of an exciton (X = e + h) bound to an additional 
electron^S The trion binding energy is 

A = E X -E X -, (1) 

where E is the X/X~ ground state energy. The lowest 
LL energy £o is set to zero for each particle. 

The exhaustive reviews of experimental and theoretical 
work on trions have been written by Peeters et al., 7 Bar- 
JosephjS and also by us£ Let us hence only summarize 
here that despite involving only three particles, quan- 
tum dynamics of a trion in the presence of confinement 



and magnetic field is neither trivial nor completely un- 
derstood. The complications include the competition of 
several (cyclotron, Coulomb, subband, and Zeeman) en- 
ergy scales, complex structure of the valence band, spin- 
orbit effects (and their effect on relaxation), and coupling 
to free carriers, lattice defects, and phonons. 

However, it is established that, depending on the pa- 
rameters (composition and width w of the quantum well, 
electron concentration g, magnetic field B, etc.), the trion 
energy spectrum contains one or more bound states, dis- 
tinguished by the total spin of the pair of electrons S 
and the relative angular momentum M.. (i) At small B, 
the only bound trion is the "singlet" X~ with 5 = 
and M. — 0, a 2D analog of the Hydrogen ion^ifi It was 
first observed in CdTe by Kheng et al& and then also in 
GaAa 11 . 12 ! 13 ! 14 and ZnSe±£J£ (with considerably differ- 
ent effective Rydbergs Ry and Bohr radii ae). Positive 
singlet trions X + = 2h + e (with the opposite effective- 
mass ratio) were also observed 1718 in acceptor-doped 
samples, (ii) In the (unrealistic) limit of very high B 
and vanishing w, the singlet trion unbinds as predicted 
from the "hidden symmetry"i2*22i2i (particle-hole sym- 
metry between conduction electrons and valence holes 
in the lowest LL). However, it is replaced by a differ- 
ent bound state: the "dark triplet" X^ with S — 1 and 
A4 = — l)22i and with an infinite optical lifetimeiSiSi De- 
spite initial difficulties 2 ^ it was also eventually detected 
in PL^SSiSLSS and its vanishing oscillator strength was 
confirmed directly in optical absorption^ The singlet- 
triplet crossover occurs in rather high magnetic fields^ 
and it is fairly sensitive to the parametersi2ii22i22i (iii) Ad- 
ditional, less strongly bound trions occur at intermediate 
fields, including the "bright triplet" X^ with S — 1 and 
A4 = 0. 31 All three trions were later confirmed by several 
independent calculations^ 2 *^ and experimentsi2L2£if2ii 

Trions were also identified beyond the "dilute regime," 
in PL spectra of the fractional quantum Hall states. It 
was founds that a trion immersed in a Laughlin v = 1/3 
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liquid may survive in form of a fractionally charged 
"quasiexciton" state which retains the characteristic in- 
ternal e—h correlations of its trion parent. This extends 
the variety of systems in which trions occur and also al- 
lows to use their recombination spectrum as an indirect 
probe of incompressible electron fluidsi^&Si This possibil- 
ity seems especially attractive for systems in which these 
fluids may form but cannot be studied by transport. 40 
Remarkably, the recombination spectra of the quasiex- 
citons formed from different trions and interacting with 
fractionally charged quasiparticles of the Laughlin liquid 
behave differently when the filling factor passes through 
v = 1/3: they show a discontinuity for the X^j but no 
anomaly for the X~. Therefore, the quantum wells for 
studying incompressible fluids by PL need be designed so 
as to yield the XT, ground state. This could be greatly 
helped by accurate theoretical predictions. 

In this paper we address this problem by extensive 
realistic numerical calculations of the binding energies 
of (negative and positive) trions in doped GaAs quan- 
tum wells in high magnetic fields. Compared to the ear- 
lier studies asymmetry of the quantum well 
caused by one-sided doping is fully taken into account 
by self-consistent calculation of the subband wave func- 
tions and by inclusion of up to three subbands along with 
several Landau levels in the Hilbert space. Especially 
for the the singlet trion and for doped wells, the results 
are significantly different from the previous, less accurate 
models. 

Another goal of this paper is to explain and justify the 
configuration-interaction exact numerical diagonalization 
in Haldane spherical geometry for the realistic calcula- 
tions of the few-body excitonic complexes. The method 
was recently used in several related studies without giv- 
ing a detailed descriptionj^Lli leaving the questions of 
its validity and accuracy largely unanswered. The de- 
tailed Sees, [n] and IIHI fill this gap, and with the future 
studies in perspective, point out the advantages over the 
(complementary) variational approch. 



II. MODEL 
A. Haldane spherical geometry 

In Haldane's spherical geometry^ the particles (elec- 
trons or holes) are confined to the surface of a sphere of 
radius R. The radial (i.e., normal to the surface) mag- 
netic field B is produced by a Dirac monopole placed at 
its center. The monopole strength 2Q = AttR 2 B/c/)q is 
defined in the units of the <po — hc/e quantum as the 
total flux through the surface. This definition can also 
be rewritten as R 2 = QX 2 which makes clear that Q 
measures the surface curvature in the units of magnetic 
length A = y/ hc/eB (rather than defining a specific B). 
Also, due to Dirac's condition, 2Q must be integral. 

The single-particle states iQ,i,m form degenerate shells 
of the eigenstates of angular momentum length / and pro- 



jection m. They are called "monopole harmonics."— The 
nth shell has I = Q + n and it corresponds to the nth 
Landau level (LL). Its degeneracy is g = 21 + 1 (m going 
from —I to I), which scales linearly with B and R 2 . 

In the lowest LL (n = 0), the m — ±Q states describe 
the closest orbits around the opposite poles, while, in 
general, (z) oc m. The wave functions scale with A. In 
the X/R — > limit, the sequence of orbitals with m = Q, 
Q — 1, Q — 2, evolve continuously into the series 
of states with angular momenta /i = 0, — 1, —2, 
obtained in the symmetric gauge on a plane. Therefore, 
also the spherical orbitals with m = Q+fi can be denoted 
by the value of /i (both in the lowest the excited LLs). 

The spherical geometry was originally proposed for 
numerical studies of extended systems (such as incom- 
pressible quantum liquids). Like confinement or periodic 
boundary conditions, it allows for having a finite number 
of electrons (as required for the exact numerical diago- 
nalization of the hamiltonian) in a finite area (to achieve 
a finite density of the actual, infinite system). How- 
ever, Haldane's geometry does it without breaking the 
2D translational symmetry, essential for various proper- 
ties of the liquid. The 2D rotational symmetry of the 
sphere replaces the symmetry of the plane and leads to 
the conservation of (the same number of) two orbital 
quantum numbers - two components of the total angular 
momentum: length L and projection L z . 

Due to the relation between orbital Hilbert eigenspaces 
in both geometries, the "spherical" L and L z correspond 
to the "planar" quantum numbers. The neutral exci- 
tations, which on a plane move along straight lines with 
well defined wave vectors k. on a sphere move along great 
circles and carry a conserved L = kR (the degeneracy 
associated with L z corresponds to that of the orienta- 
tion of the wave vector, and the quantization of L and 
L z is a finite size effect). The charged excitations move 
along the (quantized) cyclotron orbits in both geome- 
tries, and L must be interpreted in terms of the planar 
angular momentum M rather than in terms of k (in this 
case, finite-size affects the orbit by the surface curvature 
only). Clearly for the spherical geometry to be useful, 
the characteristic length of correlations responsible for 
the studied effect must scale with A, so that it can be 
made smaller than R (and the finite size effects can be 
eliminated by the 2Q — > oo extrapolation). 

Haldane's geometry has also become useful in study- 
ing excitons or trions. The sole advantage is the 2D rota- 
tional symmetry of the sphere, preserving the degeneracy 
of the trion LLs in a finite size calculation. The decou- 
pling of trion's relative dynamics from the cyclotron mo- 
tion of its "center of mass" (CM) is a big help in the 
identification of different trions in the energy spectra, es- 
pecially of the less strongly bound states. We put CM in 
parentheses here, since it is not exactly the CM motion 
that decouples for charged complexes in a magnetic field. 

Different trion states are distinguished by their relative 
(spin and orbital) correlations, reflected in difference in 
the relative angular momentum M. . The CM excitations 
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carry an arbitrary (negative) angular momentum M cul 
and lead to each trion having its own LL accommodat- 
ing different total angular momenta M = M. + A4 cm . Let 
us stress that the proper definition of M. and M. cm inde- 
pendently conserved in the magnetic field is non-trivial 
by involving the "magnetic translation" generators'^ 

If, in order to facilitate numerical diagonalization, the 
single-particle basis is restricted to several lowest single- 
particle orbitals: \\x\ < fJ. max , then A4 cm and M. are no 
longer independently conserved. This causes artificial de- 
pendence of computed trion energies on M and makes 
distinction of different trions more difficult. This prob- 
lem does not occur in spherical geometry, in which the 
trion LLs have a form of angular momentum multiplets 
with L = Q + M. . Resolving L (a conserved number re- 
gardless of finite LL degeneracy) in the diagonalization 
makes the identification simple: each bound L-multiplct 
is a genuine trion state. Exact diagonalization in cither 
geometry may be combined with the size extrapolation to 
eliminate finite-size (qualitative and quantitative) errors. 

Let us summarize the differences of both approaches. 
On the plane, the Coulomb matrix elements are accurate 
at each value of /i ma x, but the trion LL degeneracies are 
only restored in the /z max — > oo limit. On the sphere, the 
2D symmetry holds at each 2Q, but the interaction ma- 
trix elements depend on surface curvature. The last effect 
can be reduced by calculation of matrix elements from 
Haldane pseudopotentials for a plane (using the spheri- 
cal expansion and Clebsch-Gordan coefficients), helpful 
when the 2Q — > oo extrapolation is impossible. 



B. Hamiltonian matrix elements 

1 . Basis 

For diagonalization of the hamiltonian H we choose 
the configuration interaction (CI) basis, e.g., 
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for the trion. Here and S are electron and hole cre- 
ation operators, and the composite indices i, j, and k in- 
clude all relevant single-particle quantum numbers (spin 
a, well subband s, LL index n, and angular momentum 
m). Although for only two or three particles the CI ba- 
sis does not produce as sparse matrices of the (two-body) 
hamiltonian as for many-electron systems, it is still a con- 
venient choice for the evaluation of matrix elements. The 
single-particle states used in the exciton and trion basis 
include both spin projections (a =] and |), up to five 
lowest LLs (n < 4) and three lowest subbands (s < 2) 
for electrons and holes. Different values of 2Q < 30 were 
used to eliminate dependence on finite LL degeneracy. 

The pair electron spin S (for the trion) and total an- 
gular momentum L as well as their projections S z = Y] cr 
and L z — are all conserved quantum numbers. The 
CI basis is the eigenbasis of S z and L z , while the energy 



(excluding the Zeeman term) only depends on S and L. 
Therefore, only the subspace corresponding to the mini- 
mum values of S z = and L z — (or 1/2, if half-integral) 
need be diagonalized, containing all S and L states and 
thus covering the entire energy spectrum. 



2. One-body terms 



The exact cyclotron energy on the sphere is 



, I n(n+l) 

e„ = riLu r \ n H 1 — 

1 2 2Q 



(3) 



In calculation, the last term equal to n(n+l)fr /2J (here, 
,/ = m* R 2 is the moment of inertia on a sphere, and m* 
is the effective mass, not to be confused with angular 
momentum) was omitted as a finite-size artifact, and the 
realistic values of huj c for GaAs were taken from experi- 
ment. For electrons, huj c / B = 1.78 meV/T; for the heavy 
holes we used the following interpolated formula 



TlUJr 
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with the following parameter set: a = 0.45 meV, 7 = 
0.282 meV/T, ft = 275 meVnm 2 and ft = 10 meVnm 2 , 
taken after Cole et al£L and expected to work adequately 
at B > 10 T in w = 10 — 30 nm wells (in doped wells, 
the actual effective width of the hole layer is used for w). 

Although Eq. J2J describes well only the n = — * 1 
gap, we used it (for the lack of better information) to 
model scattering to higher LLs as well, and tested sen- 
sitivity to this approximation. We also neglected the 
mixing of heavy and light holes (off-diagonal terms of 
the Luttinger-Kohn hamiltonian). We expect these ap- 
proximations to affect the negative trion binding energy 
A much less than each of the exciton and trion ener- 
gies, Ex and E x ~ , separately. However, a more accurate 
treatment of the hole LL spectrum might considerably 
improve the model for the positive trions. 

The energies of higher quantum well subbands (excited 
states in the z-direction) £ s were calculated as a function 
of well width w and electron concentration in the well g. 
In doped wells (g > 0), the electric field of the remote 
ionized donors penetrates the quantum well, affecting the 
positions of occupied energy levels in the well (for our 
parameters, it is only £q) with respect to the Fermi en- 
ergy. This in turn affects g and the sheet concentration 
of ionized donors, and thus also the electric field itself. 
The Schrodinger and Poisson equations must be therefore 
solved self-consistently to give the values of £ s 4£ 

As an example, let us look at the values for the s = 1 
subband (we set £0 = 0), which for w — 10—30 nm can be 
fitted by power-law curves £\ = ^w~ K . For the electrons 
we obtain £ = 3600 and k — 1.60, weakly dependent 
on doping, and for the holes, ( — 1400 and k = 1.87 
for a symmetric well and C — 430 and k = 1.25 for 
g = 2 • 10 11 cm~ 2 (the values appropriate for £ measured 
in meV and and for w measured in nm). Clearly, strong 



subband mixing can be anticipated at least in higher 
fields and in wider wells (e.g., compare £\ = 30.5 and 
6.9 meV to %ui c = 44.5 and 7.6 meV for electron and hole, 
respectively, in a symmetric 20 nm well at B = 25 T). 

Since the projection of pair electron spin, S z , is con- 
served by the Coulomb interaction, the Zeeman energy 
Ei does not have to be included in the diagonalization. 
Instead, the entire energy spectra obtained from the diag- 
onalization without Ez must be at the end rigidly shifted 
by S z Ez- For the trion binding energies it means sub- 
tracting one electronic Ez from A of the singlet states 
(as needed to flip the spin of one of its two unpolarized 
electrons). In this approximation, the spin-orbit coupling 
is neglected and it is assumed that that electron and hole 
Ezs are both constants. In fact, they depend on n and 
s, through Rashba and (in asymmetric structures) Dres- 
selhaus couplings, both mixed by Coulomb interaction^ 

In the recombination spectrum of excitons or trions, 
the sum of electron and hole Zeeman energies splits the 
"±1" peaks corresponding to two circular polarizations 
of light (both present if both electron and hole spins 
are populated). If Ez were constant, the same splitting- 
would result for the X and all X~ states. The observed 
differenc e) 18 ! 50 reveals the inaccuracy of this approxima- 
tion. It also reflects different LL and subband mixing 
in each of these complexes. In our model, the spin-orbit 
effects can be included as the zeroth order perturbation, 
given n- and s-dependence of Ez (the LL and subband 
occupation numbers are easily extracted from the X and 
X~ wave functions obtained from the diagonalization). 



where we use standard notation for the reduced function 

< m Mi) = (-i) p E^;-^ (v«i) 2r+?) (9) 

r 

with u\ — cos(#i/2), v\ = sin(0i/2), p = f] — m. 



I + r)\ (I — rj\ (I + m\ (I — m 

T ) V T / \ T +P 

(10) 

and the summation going over all r for which none of the 
arguments of the factorials in £ are negative. 

In Eq. (J5J, we take the chord distance r = (fii ~~ ^2)R- 
From now on, let us use 9 = 9\ 2 for the angle between 
fii and il 2 . The Coulomb potential becomes 



V(6) 



2eRv 



(11) 



By applying rotation T> under the integral over fii we get 
(n'^m'^n'z, m' 2 \v^ 1 \ni,mi;n2, m 2 ) (12) 

= ^ m 'l\ V l 1 2 .*n' 2 V 'n 2 ,rn2 K> m l) K^l V ^ \ n 2,V) ■ 

n 

Using Eqs. l|Sl- llUfl . the first matrix element above is 
(n\,m\\V 1 ' 2 * ,Vl 2 m K,toi) 

_ rm'i+m'z f l' 2 ,»' f l 2 ,a r™i,mi / .\ /-, o\ 

" "mi+mj 2-^i^>ri,rn' C *ri,rn2 ln i> m i\ Uj -< u l V 10 ) 



3. Two-body terms 

a. Ideal 2D layers. Let us begin with derivation of 
the matrix elements of the Coulomb interaction V ee = V, 



V(r) = - 
er 



(5) 



(with e being the dielectric constant) for the layers with 
zero width. We use rays in spherical coordinates fi = 
(9, (f>) as the arguments of monopole harmonics and write 

(a[, n'i, rn'i, a' 2 , n' 2 , m' 2 \V\ai, ni, mi; a 2 , n 2 , m 2 ) 



s^Xi d 2 n lC i 2 n 2 Y* 



,(fi 2 ) 



xy(r)y g , il>mi (fii)y Q , i2>m2 (fi 2 ). (6) 

This four-dimensional integral is solved by using the ro- 
tation relations for monopole harmonics^ 

T>Y Qtl>m (n) = e^J2Km Y Q^(.ty, (7) 
n 

different from those for standard spherical harmonics 
only by an inconsequential phase factor ip(Q). The Euler 
rotation 45 T>(0, —0i, —<f>i) moves fii to the north pole of 
the sphere N . Its Wigner function is 

Km = Km(0, -Ox, -th.) = e <m<Al < ro (-^i), (8) 



with b = 2(a! + a + -q) ~~ m' 2 — m 2l a = 2(1' 2 + 1 2 ) — b, and 
&\a,b) = 2Trfdcos9 Y Q , v , m ,{6) u a v b Y Q , t , m (d). (14) 



The second integral in Eq. (|12fl is, up to a prefactor, an 
effective one-body Coulomb pseudopotential 



V™ (m) = (n',m\V(9)\n,m) 



(15) 



describing scattering of an electron or a hole by a point 
charge (e.g., ionized impurity) located at the north pole. 
Note that it is different from the definition of Haldane 
(pair) pseudopotential. The V(m) for the sphere can be 
converted into V(fi) for the plane via \x = m — Q. This 
interaction conserves m, but allows for an inter-LL tran- 
sition n — > n' ' . Using the symbol / it simply reads 

(t^,» 7 |i;- 1 |n a ,»j)=^(0,-l). (16) 
To calculate the integral I of Eq. (|14f> let us recall that 
Y Q ,i, m (e) = N Q>ltm u Q+m v Q - m pQ- m ' Q+m (cos0), (17) 
where 
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FIG. 1: (color online) Lowest electron and hole subband wave 
functions Xs( z ) m GaAs quantum wells of width w = 10 nm 
(a) and 20 nm (b) doped on one side to electron concentra- 
tion g = 2 ■ 10 11 cm -2 . Squared electron-hole overlaps are 
indicated. Note only two bound electron subbands in (b). 



and the Jacobi polynomial P can be expanded as 

P^{x)=Y J Cniu 2{n - q) V 2 \ (19) 



with the coefficients 



<i = (-i) 9 ( 



n + cA (n + j3 
n-q)\ q 



(20) 



and the summation range defined by the above factorials. 
Finally substitution of Eqs. (|17H20|) into Eq. {H} gives 



x 



(21) 



E Q-mt^+ni! Q-m 2 ,Q+m 2 at—: 



'"2,92 
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b+ (rii + n 2 )/2 and rh 

qi-q 2 )/2. 



- b+ (mi 



where Q = Q- 
m 2 +ni+n 2 

By collecting all above equations, a closed expression 
for the e-e matrix element in a layer of zero width is 
finally obtained. The e—h elements are found by particle- 
hole conjugation (with i below being composite indices) 



(ii;i 2 \ v eh\k;i2) = - (i' 1 ;i 2 \V ee \ii;i 2 ) 



(22) 



b. Layers of finite width. In the ideal 2D limit, the 
interaction matrix V depended on B only through the 
overall scale e 2 /A. A finite width w of the quantum well, 
causing finite extent of wave functions in the normal (z) 
direction, introduces an additional parameter w/ A oc \[B 
that affects matrix elements individually. This is true 
even in the simplest case when the quantum well asym- 
metry and the difference between electron and hole sub- 
band wave functions can be neglected. Hence, the matrix 
elements no longer explicitly depend only on 2Q (surface 
curvature), but also on the form of subband wave func- 
tions and on B (via the w/X length scale ratio). 

The realistic electron and hole subband wave functions 
Xs{z) for asymmetrically doped quantum wells were cal- 
culated self-consistently^ (together with the subband en- 
ergies £ s of Sec. Ill B"2|) . An example is shown in Fig. 



FIG. 2: (color online) Examples of one-body pseudopoten- 
tials (interaction energy V vs angular momentum /i) defined 
by Eq. 12411 and used in calculation of two-body matrix el- 
ements for e-e and e-h Coulomb scattering in/between dif- 
ferent LLs n (a) and subbands s (b), calculated for one-sided 
doped 20 nm GaAs quantum well at magnetic field B = 15 T. 



for two GaAs quantum wells of different widths, both 
doped on one side to the same electron concentration 
g = 2 ■ 10 11 cm~ 2 . The charge separation due to an elec- 
tric field inside the well is significant for the w = 20 nm 
well and negligible for w = 10 nm. This is also confirmed 
by the indicated squared e-h overlaps |(x e |x/i)| 2 - 

To include the finite quantum well width in Coulomb 
matrix elements, the bare potential V(r) must be re- 
placed by an effective in-plane potential containing the 
appropriate form factors, 



dz\dz 2 - 



Xs> ( Z l) Xl> (^) Xst (Zl) Xs 2 {Zl) 



yV 2 + (zi - z 2 ) 2 



(23) 

On a sphere, Vs^s 2 2 {&) is obtained from the above by re- 
placing r with 2Rv, exactly as in Eq. 1)11(1 for the ideal 
2D layers. In the matrix element (i[;i' 2 \V\ii;i 2 ) the first 
integral in Eq. (|12|l in unchanged, and the modified po- 
tential V(ff) enters only the one-body pseudopotential V, 
which is hence no longer given by Eq. (|16f) . However, 
since Xe 7^ Xh, different effective potentials V and pseu- 
dopotentials V must be used for e-e and e-h interactions. 

Not being able to calculate V analytically, we are forced 
to compute the entire tables of these 3D integrals for each 
choice of subband wave functions x an d 2Q, 



{ri,m\V££(0)\n,m) 



(24) 



The above pseudopotential describes inter-LL scattering 
of an electron or a hole by an artificial charge that is 
localized at the north pole within the 2D sphere and un- 
dergoes transition between the subband states si and s[ 
along the normal direction. It also can be converted into 
V(fi) for the plane via fx = m — Q. 

Several examples of V(/i) are plotted in Fig.[3to illus- 
trate the dependence on subband and LL indices. Despite 
a significant charge separation apparent in Fig. ^a), the 
corresponding e-e and e-h pseudopotentials are rather 
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FIG. 3: (color online) Pseudopotentials similar to Fig.|2] but 
defined by Eq. 1251 and describing interaction of an electron 
e~ or a hole h + with a point charge p located in the center 
(a) or at either wall (b) of the quantum well. 

similar. However, the difference reaching almost 1 meV 
(for w = 20 nm, g = 2- 10 11 cm" 2 , and B = 15 T) is 
sufficient to have a big effect on the trion binding energy. 

Scattering by a point charge (ionized impurity) at a 
distance d from the middle of the well is described by a 
different pseudopotential, 

KnU m ) = ( n 'MV s s : d (8)\nM, (25) 
with the following effective potential, 

Vl'Ar) = ~ I dz XM*>*«(*) ) (26) 
' e J ^r 2 + (z- d) 2 

where, as before, r is replaced by 2Rv for the sphere, and 
conversion to V(fJ>) is immediate. 

The examples in Fig. [31 illustrate the sensitivity to the 
impurity position d and to the electron or hole orbital 
indices s and n, as well as the electron/hole asymmetry. 
Remarkably, the energy scale of carrier-impurity inter- 
action is not much (though nevertheless by up to 2 meV) 
greater than of the carrier-carrier interaction of Fig. [2] 

The computation of a great number of needed integrals 
(|24|l can be accelerated by an adiabatic approximation in 
which the 3D integration is divided into the averaging of 
a squared relative z-coordinate 

{Di\\if = (4,41(21 - z 2 ) 2 \ Sl ,s 2 ) , (27) 
and integration (only over 9) of an in-plane potential 

V££(r)= , ^ , ■ (28) 

vW(££,1 2 2 ) 2 

As its previous versions ^L^L the adiabatic approximation 
relies on D 2 being small compared to (n! , m|r 2 |n, m), so 
its accuracy depends on the in-plane quantum numbers. 
We found that, for a given combination of s's and n's, 
it is satisfactory for sufficiently small m. However, the 



FIG. 4: (color online) Effective e-e, e—h, and h—h separations 
D in direction z normal to the quantum well plane defined 
by Eq. 12911 as a function of relative angular momentum /i. 
Horizontal lines mark averages D of Eq. 1271 . Different LL 
combinations (nn 1 ) are shown, for the lowest subband only, 
(a) Same system as in Fig.|2] (b) symmetric or undoped 10 nm 
quantum well. Magnetic field B = 15 T in both frames. 



exact 3D integration must still be carried out for those 
several largest m's (e.g., we chose Q — m = —fi < 5). 

A correct effective average separation D in the normal 
direction (such as to give the exact value of V) can be 
denned for each set of s's, n's, and m by the requirement 
that (for this particular choice of quantum numbers) 

Vt'llttn (m) = (n',m\ Ltk = = \n,m) . (29) 

/ 2 i ( n s i ,s 2 > n \2 
Y ~r v-L^si ,S2 ; n,m ) 

Certainly, D s s lf s llZ,m ~ -Dsl'S for sufficiently small m, 
while it begins to depend on n' , n, and m as m increases 
and approaches Q (i.e., as u = m — Q approaches zero). 
The typical dependences are shown in Fig. with the 
averages of Eq. (|27|) also marked by horizontal lines. Ev- 
idently, D depends on n (or s, not shown) much less than 
V does. On the other hand, using a constant value D for 
all values of /i is not a justified approximation. 

When the whole table of 3D integrals V is known, the 
corresponding table of D's can be immediately calculated 
by numerical solution of the ID equations l|29|l . We found 
that the values of -D(u) do not depend on the surface 
curvature nearly as much as V(/i). Therefore, for a given 
quantum well (the width w and the set of subband wave 
functions x) an( i magnetic field B, the time-consuming 
exact 3D integration of the V table need only be per- 
formed for one value of 2Q and for several largest values 
of u (e.g., we used 2Q = 20 and u > —5). The approx- 
imate (but accurate) tables of V for other values of 2Q 
are then quickly recovered from the corresponding table 
D, quickly computed from either Eq. (|27|l or 1)29(1 . de- 
pending on fi. Finally, the two-body Coulomb matrix 
is calculated (analytically) from the coefficients V in the 
same way as shown in Sec. Ill B 3 bl for the ideal 2D layer. 
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C. Exact diagonalization 

The hamiltonian matrices H were diagonalized numer- 
ically using a Lanczos based algorithm. The matrix was 
stored in a packed "compressed-row" format and divided 
into a number of segments that could fit in computer's 
memory. The segments were loaded one by one during 
the matrix-vector multiplication (the only operation per- 
formed on the matrix - once in each Lanczos iteration). 
Since the vanishing of a matrix element (i',j';k'\H\i,j; k) 
depends only on the two basis states, the location of non- 
zero elements (the most time-consuming part) had to be 
only done once for all used combinations of w, g, or B. 

To resolve all quantum numbers conserved by H and 
obtain the energy spectrum in the form of E{S Z , S,L Z ,L), 
the Lanczos procedure must be carried out inside the 
eigensubspaces of pair electron spin S (for the trion) and 
total angular momentum L (note that S z and L z are 
automatically resolved in the CI basis). The matrix ele- 
ments of ladder operators S + = J2 a+ an d L + = ^2,l + 
(summations go over all particles) have a simple form 
in the CI basis, allowing analytic expression of the L 2 = 
L+L-+L Z (L Z -1) and S 2 = S+ S~ +S Z (S Z -1) matrices. 

The initial Lanczos L«S-eigenstates with desired eigen- 
values of L 2 and S 2 were calculated using a conjugate- 
gradient method^ To avoid dealing with two indepen- 
dent eigenproblems, we used a joint operator S 2 + 2L 2 , 
whose eigenstates are simultaneous and unique eigen- 
states of both L and S (owing to their discrete spectra). 

The accumulation of numerical errors leading to escape 
from the initially set LS subspace after more than a few 
Lanczos iterations was remedied by the LS'-projection of 
each consecutive Lanczos state. When more than one 
lowest eigenstate was needed at some L and S, orthog- 
onalization to all preceding Lanczos states was done at 
every iteration to avoid "ghost states." Only three most 
recent Lanczos states were held in computer's memory at 
the same time; all earlier ones were stored on a disk, in 
need of either the above-mentioned orthogonalization or 
recalculation of the eigenstates into the original CI basis. 

Depending on spin and orbital symmetry, the diago- 
nalization need not be repeated in every [S z ,S,L z ,li) 
subspace. In the absence of an impurity, it is only done 
for L z = (or 1/2 if half- integral) . Furthermore, the 
spin quantum numbers are only relevant for the trions, 
and their calculation was only carried out for S z = 0. 



III. CONVERGENCE AND ACCURACY 

Being essentially a perturbativc method, the CI exact 
diagonalization converges as a function of the size of the 
many-body CI matrix. Fast convergence requires that in- 
teractions are weak compared to the single- particle gaps. 
On the other hand, the maximum size of the CI matrix 
is limited by both the knowledge of low-energy single- 
particle spectra and by the computational capability. 



TABLE I: Dimensions of the 2e + h configuration- interaction 
bases and corresponding numbers of above-diagonal non-zero 
matrix elements of the hamiltonian H, squared angular mo- 
mentum L 2 , and squared 2e spin S 2 , for the S z = L z = 
subspace and for different monopole strengths 2Q and maxi- 
mum subband and LL indices s max and n ma x- 



2Q 


Smax 


"•max 


dimension 




N H 


N L 




N s 


20 








331 




7620 


2191 




640 


20 





4 


0.6 • 10 5 


0.4 


■ 10 s 


0.4 ■ 10 6 


1.1 


10 5 


30 





4 


1.1 ■ 10 5 


1.1 


• 10 s 


0.8 • 10 6 


2.3 


10 5 


20 


1 


4 


4.6 ■ 10 5 


13.5 


■ 10 8 


3.1 ■ 10 6 


9.2 


10 5 


20 


2 


2 


2.9 • 10 5 


6.5 


■ 10 s 


1.9 ■ 10 6 


5.7 


10 5 



(a) w=20nm, p=2.5x10 11 crrr 2 , 

B=15T 



(b) w=10nm, symmetric, B=25T 
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FIG. 5: (color online) Dependence of binding energies A of 
different trions X - (distinguished by the 2e spin S and rela- 
tive angular momentum M) on squared surface curvature (R 
is the sphere radius; A is the magnetic length), calculated in- 
cluding five LLs and the lowest subband only for (a) one-sided 
doped 20 nm GaAs quantum well at magnetic field B = 15 T, 
and (b) symmetric/undoped 10 nm well at B = 25 T. 



In Tab. [I] we list sizes of the H, L 2 , and S 2 matrices 
for several CI bases used in computation. The 2e + h 
hamiltonians are not very sparse (2iVjy/dim > 1%), so 
numerical feasibility is defined almost exclusively by the 
number of their non-zero elements Nh (rather than by 
the space dimension or the sparsity of L 2 or S 2 ). For the 
largest basis we have used, compressed storage of a single 
i/-matrix required over 15 GB of disk space. 

Constrained by the manageable-size restriction and by 
limited knowledge of the matrix elements, we are still able 
to calculate some quantities with useful accuracy and to 
model various dependences at least qualitatively. Below 
we present several tests carried out to estimate conver- 
gence of our calculations of the trion binding energies A. 

First, we analyze the finite-size and surface-curvature 
errors. In Fig. [5] we plot A calculated in two different 
systems (a narrow symmetric well in a high magnetic 
field and a wider asymmetric well in a lower magnetic 
field), including one subband and five LLs, as a function 
of the surface curvature (A/i?) 2 = Q . The dependence 
is fairly weak, justifying the use of spherical geometry 
for the quantitative calculation. Moreover, the regular 
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FIG. 6: (color online) Similar to Fig. but with dependence 
on inverse maximum LL index n max (number of included LLs) 
calculated for 2Q = 20 and including the lowest subband only. 



behavior of A(<3 -1 ) allows reliable extrapolation to the 
(planar) R — > limit. In the following we used 2Q = 20. 

In Fig. we study dependence on the number of in- 
cluded LLs. Although electron and hole cyclotron en- 
ergies are different, we used the same maximum value 
n majt for both particles, because the rotational symme- 
try of a sphere (or the magnetic translational symme- 
try of a plane) forbids mixing of different electron and 
hole LLs in the optically active (k = 0) exciton ground 
state: (n' e , m' e ; n' h , m' ?l |V~|0, m e ; 0, to/j) cx . The de- 
pendence of A on n max is much stronger for the singlet 
trion than for both triplets. This is due to the stronger 
Coulomb interaction in this most compact state (the 
spin-unpolarized electrons have a large projection on the 
pair state with the minimum relative angular momentum 
/i = 0, not allowed for in a triplet state), causing a more 
efficient mixing of different LLs. In the following we in- 
cluded five LLs (n max = 4). This promises to describe 
fairly accurately both triplet states, but A for the singlet 
state will probably be noticeably underestimated (by 0.2 
and 0.4 meV in the two shown examples). 

Inclusion of a great number of LLs is more difficult 
for the holes due to the complicated structure of the va- 
lence band. Since the correct energies e„ of higher LLs 
are not known in general (for arbitrary w, g, and B), 
we largely ignore this complication by assuming equidis- 
tant LLs separated by the cyclotron energy taken from 
experiment for the lowest (n = — * 1) transition. In 
Fig. [7| we test sensitivity of the results to this assump- 
tion by studying the dependence of A on the dimension- 
less factor ( reducing the select (nth) hole's cyclotron 
gaps Ae n = e n +i — £ n compared to Tluj c of Eq. Q}. This 
test we have done for negative and positive trions, antic- 
ipating that the X + could show more dependence. The 
thin lines were obtained for all hole LLs spaced by the 
same Ae = £hui c . A more reasonable test is shown with 
the thicker lines, obtained for only the (unknown) gaps 
beyond the n = 1 LL being rescaled by (. The lack of 
strong dependence justifies the neglect of exact LL struc- 
ture for the holes in the following calculation. 

Finally, in Fig. [S] we analyze the effect of higher quan- 




0.4 0.6 „ 0.8 1.0 0.4 0.6 „ 0.8 1.0 



FIG. 7: (color online) Dependence of negative (a) and posi- 
tive (b) trion binding energies A on factor ( reducing the hole 
cyclotron gaps compared to Eq. calculated for symmet- 
ric/undoped 15 nm quantum well at magnetic field B — 10 T, 
including five LL and two subbands. Thin lines: e„ = £o+ 
nC^hwc (n > 1); thick lines and symbols: si — £o + hu c and 
£ n — £i + (n — 1) (hu c (n > 2), i.e., only higher gaps reduced. 
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FIG. 8: (color online) Comparison of magnetic field depen- 
dences of negative and positive trion binding energies A(B) 
calculated including five LLs and between one and three sub- 
bands, (a) X~ in one-sided doped 20 nm GaAs quantum well; 
(b) X + in symmetric/undoped 15 nm well. 

turn well subbands. On two examples (negative trions in 
a doped 20 nm well and positive trions in a symmetric 
15 nm well) we compare the magnetic- field dependence 
of A calculated using three LLs and between one and 
two subbands. Analogously to the LL mixing, by far the 
strongest subband mixing occurs for the singlet trion. 
The effect is especially strong in both systems used here 
as the examples due the strong charge separation in the 
lowest subband for the X~ and having two heavy-mass 
particles in the X+. Evidently, the s = 1 subband must 
necessarily be included in the calculation for all trions in 
the presence of a doping layer, and at least for the sin- 
glet state in symmetric or undoped wells. However, the 
s = 2 and higher subbands can be neglected (except only 
for the singlet trions in wide asymmetric wells in high 
magnetic fields). Hence, two subbands (s max = 1) were 
included in the following calculations. 

The results of this section can be summarized by two 
conclusions that we believe to be important in the com- 
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parison of numerical and experimental studies of trions. 
First, the dark and bright triplet states can be accurately 
modeled by exact diagonalization. The errors due to fi- 
nite size, limited CI basis, or uncertainty in some micro- 
scopic parameters either do not affect these trions or can 
be eliminated (without need for extraordinary computer 
power). Second, numerical calculations for the singlet 
trion (by the CI exact diagonalization) are far more diffi- 
cult due to a stronger sensitivity to the unknown sample- 
dependent parameters. However, note that even for the 
singlet X + we reached quite satisfactory agreement with 
the experiment of Vanhoucke et al. (cf. Fig. 3 of Ref.l28ch 

A. Variational treatment of higher subbands 

As a complementary account of the subband mixing 
we also used the following combination of exact diago- 
nalization (in the plane of the well) and the variational 
approach (in the normal direction). The pair of model 
lowest-subband electron and hole wave functions 

X(z) = ^2/w*cos(zn/w*), (30) 

where w* = w + 3.3 nm and — w + 1.75 nm denote the 
effective layer widths adequate for GaAs/Aln.35Garj.65As 
symmetric wells and for w = 10 — 30 nm, were each 
rescaled by an arbitrary shrinking factor £ < 1 

xdz)=C 1/2 x(C 1 z)- (31) 

The diagonalization was then carried out separately for 
the e + h and 2e + h systems and for different com- 
binations of electron and hole trial subband states X£- 
The result are the energy maps E(^ e ,^h), separate for 
the exciton and each trion state. Having ^ ^ 1 costs 
single-particle energy, but stronger confinement imposed 
by £ < 1 decreases total Coulomb energy of each bound 
state. The interpolation and minimization of E with re- 
spect to both variational parameters £ yield, indepen- 
dently for the X and each X - , the ground-state sub- 
band parameters (£ e ,£/i) arL d energies f?(£ e ,£h). The 
variational trion binding energies were calculated from 
Eq. (JIJ and, finally, the correction 6 A beyond the lowest- 
subband approximation was obtained by their compari- 
son with the reference values obtained for £ e = ^ = 1. 

In two frames of Fig. [5| we plot the calculated £ e and £/j 
as a function of w and B. The inaccuracy of the lowest- 
subband approximation, manifested by the interaction- 
induced distortion (squeeze) of the subband wave func- 
tions and measured by £, clearly increases as a function 
of both B and w. Both £ e and £/, are roughly linear in 
(w/\) 2 ex w 2 B. However, £/, is considerably smaller than 
£ e > 0.98 (which further enhances the < w* asym- 
metry). Both £'s differ between the exciton and different 
trion states. Remarkably, £/, and £ e are anti-correlated 
in a sense that £/j decreases and £ e increases when going 
from X to X^, X^, and X~ (at any given w and B). 

The squeeze of subband wave functions shown in Fig.[!5] 
confirms the main conclusions drawn from inclusion of 
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FIG. 9: (color online) Variational coefficients £ of Eq. 131L 
rescaling widths of electron and hole wave functions compared 
to the noninteracting lowest-subband eigenstates, calculated 
for the exciton and trions in an undoped/symmetric GaAs 
well as a function of magnetic field B (a) and width w (b). 
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FIG. 10: (color online) Variational correction SA to the trion 
binding energies due to the squeeze of electron and hole wave 
functions compared to the noninteracting lowest-subband 
eigenstates, calculated for undoped/symmetric GaAs wells as 
a function of magnetic field B (a) and the well width w (b). 

higher subbands directly into exact diagonalization (cf. 
Fig. |SJ). It predicts dependence on w and B, and ex- 
plains why the subband mixing affects the singlet trion 
much more than both triplets. The variational correction 
SA beyond the lowest subband approximation are plot- 
ted in Fig. Clearly, the quantitative agreement with 
exact diagonalization is rather poor (even if the compar- 
ison is made for the same n max = 2 restriction in each 
method). The reason is probably the combination of (i) 
the neglect of correlations in the normal direction within 
the variational calculation and (ii) the rather unrealistic 
choice of the trial subband wave function xz ■ 



IV. RESULTS AND DISCUSSION 

A. Exciton energy dispersion 

Let us begin the discussion of trion binding energies A 
with an obvious observation that they strongly depend 
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FIG. 11: (color online) Comparison of excitonic dispersions 
(energy E vs wave vector k) obtained with and without higher 
LLs and subbands, in different GaAs quantum wells. Insets: 
exciton Coulomb mass relative to the ideal value of Eq. (1351 . 



on the energy spectrum of the exciton being captured by 
a free electron. This spectrum consists of a continuous 
in-plane dispersion E(k) in each quantized subband. 

In an ideal 2D system (lowest LL and zero width), 
exciton dispersion is known exactly^S^i 



(32) 



with k — k\/2 and Iq being the modified Bessel function 
of the first kind. For small k, 



E(k) 



A 



1 



(33) 



The curvature of E{k) can be attributed to a "Coulomb 
mass" mx(k) = K 2 (d 2 E / <9fc 2 ) -1 , which must not be con- 
fused with the single-particle electron or hole cyclotron 
masses. For an ideal 2D system, it follows from Eq. 132|) 



mx{k) = h 2 



2 



1 



9 



(34) 



At zero wave vector, the relation between mx = mx(0) 
and the interaction strength is particularly simple, 



2h z 



m x = 



X 2 \E(0)\ 



(35) 



The magneto-exciton spectra in real quasi-2D struc- 
tures were studied by Lozovik et al.^ so here we concen- 
trate on the finite- width and LL/subband-mixing effects, 
illustrated in Fig. ^2 for two different quantum wells at 
different magnetic fields. The ground state energy E(0) 
which defines the maximum e—h attraction noticeably de- 
pends on both LL and (in the asymmetric well) subband 
mixing. In the insets we plot mx divided by the value 
given by Eq. I|35|l for an ideal 2D system. In both quan- 
tum wells we find significant enhancement of mx due to 
finite well width (by a factor ~ 2). Inclusion of higher 
LLs causes additional increase by another factor of ~ 2, 
while the subband mixing has a much weaker effect. 
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FIG. 12: (color online) Dependence of the trion binding ener- 
gies A on magnetic field B, for undoped or symmetric GaAs 
quantum wells of different widths w. Different bound trion 
states are distinguished by the 2e spin S and relative angu- 
lar momentum M. Five LLs and up to two subbands were 
included in exact diagonalization on a sphere. 



B. Trion binding energies 

Let us now turn to the presentation of our main numer- 
ical results, that is to the dependence of the trion bind- 
ing energies A in one-sided doped GaAs/Alo.3sGao.65As 
quantum wells on w, g, and B. The g — data are 
also adequate for the symmetrically doped wells. The 
calculations were carried out for 2Q — 20 and including 
two subbands (s max = 1) and five LLs (n max = 4). The 
accuracy of this basis was discussed in Sec. IIHI 

Plotted values of A are the Coulomb binding energies. 
The appropriate Zeeman terms must be additionally in- 
cluded to determine the absolute trion ground state or 
the energy splittings in the photoluminescence spectra. 

Depending on the parameters, we found up to four dif- 
ferent trions with significant binding energy in the energy 
spectra similar to Fig. 1 of Ref. |3lJ. Besides the "sin- 
glet" X- with (S,M) = (0,0), "dark triplet" X w (1,-1), 
"bright triplet" X^ (1,0), an additional "dark singlet" 
X^j (0,-2) was identified in some systems. Their binding 
energies will be denoted by A s , A td , A t b, and A sc j. 

In the following figures, the bright and dark states are 
distinguished by solid and dashed curves. On the other 
hand, the singlets and triplets are marked by open and 
full dots, respectively. The sets of curves obtained in the 
lowest-subband approximation (LSA) are also drawn for 
comparison with thin lines without the dots. 
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FIG. 13: (color online) Similar to Fig. CU but for one-sided 
doped w = 10 and 20 nm GaAs quantum wells of different 
electron concentration g. 



1. Dependence on magnetic field 
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FIG. 14: (color online) Similar to Fig. 1121 but with trion bind- 
ing energies A plotted as a function of electron concentration 
g, for a 20 nm GaAs well at fields B = 10 T (a) and 40 T (b). 
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We begin with the magnetic-field dependence shown 
in Fig. 1121 for four symmetric quantum wells of differ- 
ent width and in Fig. ^5] for two doped wells of different 
width and electron concentration. The singlet remains 
the most strongly bound state in all frames, with the 
singlet-triplet crossing pushed to much higher magnetic 
fields than anticipated by the earlier LSA calculations. 
Especially in the wide (w = 20 nm) doped wells the LSA 
fails completely, with the actual A S (B) far more resem- 
bling that in the symmetric well. In all systems, A s (£>) 
and Atb(-B) rise at weak magnetic fields, but the increase 
saturates at B = 15 — 30 T (depending on w and g), be- 
yond which both A S (B) and Atb(-B) remain essentially 
flat at A s = 1.5 - 2.4 meV and A tb = 0.2 - 0.4 meV 
depending on w and g. In contrast, A t d(£?) largely re- 
tains the characteristic \f~B dependence of an ideal 2D 
system (w = and lowest LL). The weakly bound dark 
singlet emerges in wider wells and at higher fields, with 
a monotonically increasing A s <i(-B). 

2. Dependence on electron concentration 



FIG. 15: (color online) Similar to Fig. 1121 but with trion bind- 
ing energies A plotted as a function of width w of a symmetric 
GaAs quantum well, at fields B = 15 T (a) and 25 T (b). 



dence on g. However, it must be realized that if A were 
to be understood as the binding energy of a well defined 
bound state, at most weakly perturbed by the surround- 
ing electrons, then B and g must be restricted to the 
values describing a sufficiently "dilute" system. A con- 
venient measure is the LL filling factor v — 2irg\ 2 (pro- 
portional e.g. to the area of the dark triplet trion divided 
by the average area per electron). While for B = 40 T 
the system can be regarded dilute v < 1/3 throughout 
the frame, at B = 10 T the lowest LL fills completely 
(u = 1) already at g = 2.4 • 10 11 cm" 2 . 

It is clear from Fig. ^] that the LSA fails completely 
to describe the effect of charge separation in asymmetric 
wells. It appears that, except for the extreme cases of 
highly doped wide wells, the effect can be simply ignored. 
However, if it need be included, the calculation must take 
into account the mixing with at least one higher subband. 



In Fig. E| we study in more detail the dependence on 
electron concentration g. It was clear from Fig. 1121 that 
the effect is negligible for narrow wells, so we only plot 
the data for w = 20 nm (with a relatively low B — 10 T 
and a very high B = 40 T compared in the two frames). 
Of all states, only the singlet shows a noticeable depen- 



3. Dependence on quantum well width 

In Fig. E] we plot A as a function of the width w (of 
a symmetric well) at two different magnetic fields. The 
LSA works remarkably well for the two triplet states, 
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FIG. 16: (color online) Same as Fig. 1151 but for positive trions. 



neither of which shows a significant A(w) dependence 
throughout the wide range of w = 10 — 30 nm. In con- 
trast, and in good agreement with the experiment of Van- 
houcke et al. (cf. Fig. 5 of Ref. 128k). A s (w) decreases 
markedly in the same range, although nearly not as much 
as predicted by the LSA. 

In Fig. E| we use the width dependence A(w) to il- 
lustrate the difference between negative and positive tri- 
ons. Remarkably, at B > 15 T dark triplet is the most 
strongly bound X + throughout the entire shown width 
range w — 10 — 30 nm. Only the singlet positive trion 
has a lower binding energy than its negative counterpart. 
All other trions enhance their binding when going from 
X~ to X + . This is especially true of the dark singlet 
which emerges as an additional robust bound state in 
the X + spectrum. Also, the width dependence for the 
pair of triplets is more pronounced in the X + . Similarly 
as in X - , the LSA works well for all positive trions except 
for the bright singlet. 



4- Singlet-triplet crossing 

From the comparison of binding energies of differ- 
ent trions one can try to establish the absolute trion 
ground state in a given quantum well (w and g) at a 
given field B. This might be especially important in 
the context of "quasiexcitons" ?2£ proposed to form in in- 
compressible electron fluids, such as the Laughlin frac- 
tional quantum Hall state at LL filling factor v = 1/3. 
Only those quasiexcitons formed from a dark triplet trion 
were found to cause discontinuity of the PL emission en- 
ergy at v — 1/3, suggesting that an appropriate struc- 
ture design (to achieve significant occupation of the dark 
triplet state) might be essential for optical detection of 
the Laughlin incompressibility. 

In Fig. El we plot the critical magnetic field B c at 
which the singlet-triplet crossing occurs in the trion 
Coulomb spectrum, as a function of both w and g. An 
accurate estimate of B c is quite difficult due to accumu- 
lation of errors in the two compared A's. Nonetheless, 
it is clear that inclusion of subband mixing favors the 



FIG. 17: (color online) Critical magnetic field B c for the cross- 
ing of Coulomb binding energies of the singlet and dark-triplet 
states of the negative trion, as a function of electron concen- 
tration n (a) and quantum well width w (b). Open circles in 
(b) mark B c for a positive trion in a w = 15 nm well. 

singlet trion, pushing B c to much higher values. 

The fact that the LSA underestimates B c was found by 
Whittaker and Shields^ (for symmetric wells). It can be 
explained by the same perturbative argument we invoked 
earlier to explain different effect of the LL mixing on the 
singlet and triplet states: by noting that Coulomb in- 
teractions are stronger in the singlet state, causing more 
efficient mixing of excited single-particle states. However, 
in contrast to Whittaker and Shields, and in agreement 
with exp eriment of Vanhoucke et al. (cf. Figs. 5 and 6 of 
Ref. l28b). we find that B c decreases as a function of w 
(e.g., leading to a crossing below 50 T in a 30 nm well). 

Other conclusions from Fig.El are the following. Dop- 
ing of even a fairly wide, 20 nm well has small effect on B c 
up to g ~ 2 • 10 11 cm" 2 (also, its decrease at higher con- 
centrations is most certainly an artifact of the s < 1 ap- 
proximation; cf. Fig. [SJ) . The crossing of Coulomb bind- 
ing energies calculated here is different from the crossing 
of the ground state energies due to additional electron 
Zeeman energy that must be added for the singlet (to 
lower B c compared to our plot). It is also different from 
the crossing of recombination energies in the PL spec- 
tra, due to dependence of the hole Zeeman energy on a 
particular exciton or trion state. 

In Fig. I17f b) we also marked B c for a positive trion 
in a symmetric 15 nm well (small circle for LSA; large 
circle for s < 1). In good agreement with the experi- 
mental reports, it is much smaller than the value for X~, 
predominantly due to the opposite effective-mass ratio. 

5. Effect of nearby ionized impurities 

Being charged, trions are efficiently localized by even 
fairly remote ionized impurities. Especially for the dark 
triplet, this might appear as a possible mechanism for 
breaking translational invariance of an otherwise perfect 
quantum well, needed for the explanation of experimen- 
tally observed optical recombination of this state. How- 
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FIG. 18: (color online) Energy spectra of the 2e + h system 
in the presence of an ionized donor D + at a distance d — w 
(a) and d = w/2 (b) from the well center (energy E, mea- 
sured relative to an exciton unbound from the D°, plotted as 
a function of angular momentum on a plane M), calculated 
for a symmetric/undoped w — 20 nm GaAs quantum well at 
magnetic field B = 15 T, including three LLs and two sub- 
bands. Lines connect the D + X~ states formed from different 
trions; the 5 = and 1 ground states are marked with bigger 
symbols; and the color labels indicate squared overlaps with 
the corresponding free-trion states | (X _ |D + X~ } | 2 . 



ever, Dzyubenko et alm& showed recently that this mech- 
anism is not efficient. The D+X^ state was found stable 
only at sufficiently large distances d of the ionized donor 
D + from the center of the quantum well, at which its 
oscillator strength remains small (~ 1% of the X~). 

We repeated this calculation in our model, but chose 
a wider, w = 20 nm quantum well and a higher field 
B = 15 T, in anticipation of a stronger effect of subband 
mixing induced by an off-center (d =/= 0) point charge. 
In Fig. [E|we compare the 2e + h energy spectra in the 
presence of a D + at d — w and w/2. On the horizon- 
tal axis is the total angular momentum converted to the 
plane, M = L z — Q (neither L nor M. are conserved in 
the presence of D + , leaving L z or M as the only good or- 
bital quantum number). The energy E(M) is measured 
relative to the configuration consisting of a ground state 
exciton unbound from a neutral donor state D° = D + + e. 
The plotted energy difference is hence the negative of 



A(M) =E X + Ed{M) -E(M), 



(36) 



which is the trion binding energy in the presence of a D + , 
generalized from Eq. JQ). The singlet and triplet ground 



states (states with the lowest E) are also indicated. 

The average trion-donor distance depends on both d 
and M. At sufficiently large distances of the donor from 
the well (in our example, at d > w), the trion is captured 
without a significant distortion of its wave function. The 
lowest state for each trion occurs at M = M. (trion's 
relative angular momentum). Each trion's A is slightly 
decreased, but the dependence on d or M is weak in this 
regime. In contrast, when the donor approaches the well, 
new trions bind and mix with X~ , X^j , and X^ . The ef- 
fect disappears at large |M|, but it is also relatively weak 
at M — 0, because the additional trions have Ai ^ 0. 56 
As pointed out in Ref-IBH. for d = w/2 (donor at the wall 
of the quantum well), A = 1.2 meV of the singlet ground 
state is only slightly smaller than A = 1.55 meV of a free 
XjT, but the triplet ground state essentially unbinds. 



V. CONCLUSION 

We used exact numerical diagonalization on a sphere 
to calculate the energy spectra of negative and positive 
trions in doped GaAs quantum wells in high magnetic 
fields. The results obtained with the inclusion of up to 
three quantum well subbands were compared with a dif- 
ferent method, combining exact diagonalization of the 
in-plane dynamics with a variational calculation for the 
normal direction. The obtained dependences of the trion 
binding energies A on the quantum well width w, elec- 
tron concentration g, and magnetic field B significantly 
improve over previous theoretical estimates. In particu- 
lar, the symmetric-well and lowest-subband approxima- 
tions are both invalidated for a wide class of realistic 
systems. Presented detailed analysis of the accuracy and 
convergence establishes the exact diagonalization on a 
sphere for future quantitative studies of excitonic com- 
plexes. For example, it can be used for the analysis of 
exciton and trion wave functions (e-e and e—h correla- 
tions, oscillator strengths, LL/subband projections, etc.) 
or more detailed studies of the interaction of excitonic 
complexes with ionized impurities and/or free carriers. 
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